Loschmidt echo and the many-body orthogonahty catastrophe in a qubit-coupled 
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We investigate the many-body generalization of tlie ortliogonality catastrophe by studying the 
generalized Loschmidt echo of Luttinger liquids (LLs) after a global change of interaction. It decays 
exponentially with system size and exhibits universal behaviour: the steady state exponent after 
quenching back and forth n-times between 2 LLs (bang-bang protocol) is 2n-times bigger than 
that of the adiabatic overlap, and depends only on the initial and final LL parameters. These are 
corroborated numerically by matrix-product state based methods of the XXZ Heisenberg model. An 
experimental setup consisting of a hybrid system containing cold atoms and a flux qubit coupled to 
a Feshbach resonance is proposed to measure the Loschmidt echo using rf spectroscopy or Ramsey 
interferometry. 

PACS numbers: 71.10.Pm,67.85.-d,85.25.-j,05.70.Ln 



Introduction. The long coherence times and the pos- 
sibility to control parameters very accurately in optical 
lattices allow to simulate exciting non-equilibrium effects 
in quantum many body systems. Particular attention has 
been devoted to the evolution of quantum stated after 
quenches [l[ in which a parameter in the Hamiltonian is 
changed either suddenly or gradually. By measuring the 
Loschmidt echo (LE), it is possible to get insight into the 
dynamical properties of the quantum-many body state. 

The LE is defined as the overlap of two wave functions, 
\'^o{t)) and \'^{t)), evolved from the same initial state, 
but with different Hamiltonians, Hq and H , 



(1) 



It measures the "distance" between two quantum states 
and serves to quantify irreversibility and chaos in quan- 
tum mechanics 0-0] • Furthermore, it can be used to di- 
agnose quantum phase transitionsi , and is also an im- 
portant quantity in various fields of physics, ranging from 
nuclear magnetic resonance to quantum computation and 
information theory. For the latter, the LE is of utmost 
importance since it measures how small changes during 
a time evolution cause decoherence and are detrimental 
for quantum information processing and storage. 

While the theoretical concept behind the LE is very 
clear, experimental protocols to measure the overlap of 
two wave functions, are challenging and rare. Experi- 
mentally, the LE and the closely related fidelity has only 
been observed in systems with limited degrees of free- 
dom 0, by coupling the system of interest to a sin- 
gle qubit, whose states probe the wavef unction of the 
system at different values of some external parameters. 
Also, while the LE for local perturbations - the X-ray 



edge singularity problem - is well understood, its behav- 
ior in generic many-body systems is poorly described 0] ■ 
Studies so far mostly focused on local perturbations and 
simple refermionizable spin chain modelsjsl-fioj. 

Here we study the interaction driven LE of a genuine 
interacting one-dimensional system, a Luttinger liquid 
(LL). One-dimensional strongly correlated systems often 
form such Luttinger liquids (LL), and can also be cre- 
ated from cold atoms Recently, the LL wavefunc- 
tion has been used to evaluate the overlap of distinct LL 
ground states 12, 1^, and optimal protocols producing 
maximal overlaps with reference states have also been in- 
vestigated 1^. In this work, we evaluate the LE within 
the Luttinger liquid (LL) description for a time depen- 
dent Hamiltonian, and verify the predictions of the LL 
model numerically on the XXZ Heisenberg chain using 
matrix-product state (MPS) based methods. Further- 
more, we propose an experimental setup to measure the 
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FIG. 1. (Color online) Schematics of the experimental setup. 
The black segments on the SQUID denote the Josephson junc- 
tions, the arrows stand for the total magnetic field. By chang- 
ing the eigenstate of the flux qubit, the total flux hence the 
total magnetic fleld changes, controlling the Feshbach reso- 
nance in the cold atomic LL. 
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LE of a Luttingcr liquid, where a flux qubit coupled to a 
Feshbach resonance is used to control the interaction in 
one-dimensional cold atom gas (see Fig. [1]). 

Model. A LL can be described in terms of bosonic 
sound-like collective excitations, regardless to the statis- 
tics of the original system. The LL Hamiltonian is given 

by [Hi 



+ «]) , (2) 



where gi{q) is the initial interaction, and ujq the energy 
of bosonic excitations. In the following, we shall assume 
that H = H{t) has the same form as Eq. but with a 
time dependent coupling. 



9q 9q{t) = 9t{<l) + ^9q{t). 



(3) 



where Agq{t) = [gi{q) — gf{q)]Q{t), and gf{q) is the final 
interaction strength. For a conventional LE Q{t) = 1 
for t > 0, but here we allow any time dependence, and 
only assume that Q(0) = and Q{t) = 1 after some 
" transition time" r [16| . 

The Hamiltonian H{t) is quadratic, and can be diago- 
nalized at any instance. Its initial and final quasiparticle 
spectra are simply given by uji/f{q) = {uJq - gf/j{q)y^'^, 
and the strength of interaction in these states is conve- 
niently characterized by the dimcnsionless LL parame- 
ters, K,/f = ^ [Uq - g^/f{q)\ I [Wq + ft//('7)] ■ 

To compute the LE, it is convenient to first diagonal- 
ize Eq. ([5]) by a standard, time independent Bogoliubov 
transformation, yielding Hq ~ '^q^gi-^iiq) b'^bq up to a 
constant shift of energy. In this basis, H{t) reads 

9{qit) , 



H = Y,^iq,t)b+bq 

1^0 



[bqb-q 



(4) 



with uj{q,t) = u:^{q) - Agq{t)^f^ and g{q,t) = 
Agq{t) ^^q-j , and the dots stand for an unimportant time 
dependent energy shift. We emphasize that our setting 
is very general [l3, and equally applies for a spinless 
fermion, boson or spin system[15|. 

Analytical results. In order to calculate the LE[18|. 
we need to determine the wave function of our quenched 
LL. This can be achieved by realizing that Eq. (j4]) cou- 
ples only pairs of states with q and —q. Consequently, 
the total time-evolution operator is the product of the 
separate evolution operators Uq{t) for all g's. Finding 
Uq{t) for a given pair of modes can be accomplished by 
realizing that the operators appearing in the Hamilto- 
nian, Ko{q) = {b+bq + b_qbtq)l2, K+[q) = b+btq and 
K-{q) = bqb-q are the generators of SU(1,1) Lie alge- 
bra. Generalizing the results of squeeze operators to the 
present case, and exploiting a faithful matrix representa- 
tion of the SU(1,1) generators [l^, the time evolution 



operator Uq{t) can finally be expressed as [20| 

Uq{t) = exp[C+(g,t)A + (g)] exp [Co (g, 0/10(9)] x 

y.cyip[C-{q,t)K^{q)+iLPq{t)]. (5) 

Here (fq{t) is an unimportant phase, and the coefficients 
CQ{q,t) and C±{q,t) can be expressed in terms of the 
time dependent Bogoliubov coefficients, Uq{t) and Vq{t), 



c+{q,t)^v;it)/u;it), CUq.t) 

Co{q,t)^-2\n{u*{t)), 
with the coefficients Uq(t) and Vq{t) satisfying 



vq{t)/u;it), (6a) 
(6b) 



idt 



Uq{t) 









i^{q,t) g{q,t) 
-g{q,t) -uj{q,t) 



Uq{t) 
Vq{i) 



(7) 



with [uqifi), Wq(0)] = [1,0]. We can now express the time 
evolution of any initial state 1^*0) as |^'(i)) = t^(Ol^o), 

with U{t)=Y{^^^Uq{t). 

In the following, for simplicity, we take [^^o) as the 
ground state wavefunction of i/o, which is the vacuum for 
the b bosons. Then the time evolved, normalized wave 
function simplifies to 



(8) 



with $(t) an overa.ll phase factor. This generalizes the 
equilibrium results 12, ll3| to the non-equilibrium situa- 
tion, encoded via the time dependent Bogoliubov coeffi- 
cients. Then the time evolution under the action of 
is trivial, since ['I'o) only picks up a phase, and the LE 
takes a particularly simple form 



C{t) — exp 



-Y,\n{\uq{t)[^ 

q>0 



(9) 



This result, in combination with Eq. ([7]) determines the 
complete time dependence of the generalized LE in a LL. 
It holds true for any non-equilibrium evolution, and ex- 
presses the LE solely in terms of the number of excited 
quasiparticles in the final state. We regularize the q sums 
in Q by an exp(— q;|(7|) factor, with 1/a an ultraviolet 
cutoff [il. 

We can now use the function Q{t) to calculate the LE. 
Changing Q{t) adiabatically, one obtains from Eq. ([7]) 
\uf{t)\'^ = 1/2 + {KjKf + Kf/Ki) /4 for times t > r. 
The LE is in this case just the overlap of the ground states 
of the initial and final Hamiltonian and, in agreement 
with Refs. 



121, ll3| , reads as 
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KKf " 





-L/2iTa 



(10) 



with L being the system size. This remains valid in the 
steady state for near-adiabatic quenches, t ^ a/v with 
V being the sound velocity in the final state. 
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For Q(i > 0) = 1, i.e., a conventional "sudden quench" 
(SQ) Um LE, however, Eq. ^ yields 



(11) 



with ujf{q) = v\q\ denoting the excitation energy after 
the quench. By plugging this back to Eq. ©, we find 
for the short time limit (t <ti a/v) 



£sqW ^exp(-cL (t/t,)Va) 



(12) 



where c is a non-universal constant of order unity. The 
characteristic decay time of this expression is tc = 
4a/v\K f / Ki — Ki/Kf\, and ^ can be identified as 
the variance of energy per particle after the quench per 
particle. For intermediate times t ^ tc, the LE displays a 
non-universal transient signal (see Fig.[2lD), however, for 
very long times, t ^ a/v, the LE for SQ becomes time 
independent and universal: 



CsQit:$>a/v) 



1 (K^ 

4 U/ 



K, 



-L / TTa 



(13) 



which holds also true for fast quenches, r <^ a/v. Exci- 
tations are only produced at i = 0, which interfere with 
each other for a short amount of time, causing Eq. ([T^. 
but after this phase coherence is lost, excitations prop- 
agate independently and only their total number deter- 
mines the overlap. 

By Eq. (|T0|) we just established that for < oo 



C 



SQ 



''ad ' 



(14) 



This can be understood in terms of simple physical ar- 
guments and must hold in general, too[lO|: the adiabatic 
LE measures just the square of the overlap {G\Go) of 
the ground state of the initial and final Hamiltonians, 
which is the probability weight of the ground state of 
H in |Go) = |$o)- For a regular LE, ($(i)|$o(i)) 
($o|e'*^e~'*^" |$o)- Inserting here the complete set of 
eigenstates of H between the two propagators, the ex- 
cited states amount in oscillating terms, interfering de- 
structively, and only the ground state contribution re- 
mains, yielding asymptotically in the thermodynamic 
limit |($(t)|$o(0)l = l(G|Go)P, and thus Csqit) = 
mt)\Mt))? ^ KG|Go)|4 = Cl^ . Using the anal- 
ogy to work statistics [l8|, this is the square of the prob- 
ability to stay in the ground state after the quench. 

The exponent of the LE is further enhanced [i^] by re- 
peating the Q{t > 0): SQ to 1, holding time t, SQ to 
0, holding time t sequence n times (bang-bang protocol). 
The generalized LE is the nth power of the SQ overlap 



in Eq. ^ as Cn{2nt) = Ci.Q{t) 



rin 



in the post- 



quench steady state, therefore the exponent is enhanced 
by a factor of n. 

Numerics. We tested the analytical predictions numer- 
ically on the one-dimensional XXZ Heisenberg modelfisj. 
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FIG. 2. (Color online) a.) The exponent of the generalized 
Loschmidt echo is shown for SQ n = 1 (circle), double quench 
n — 2 (triangle) and adiabatic time evolution (square) for the 
XXZ Heisenberg model in the steady state, starting from the 
XX point and ending up at a finite Jz, obtained numerically 
from MPS based methods (when convergence was reached). 
The solid lines are the analytical results from Eqs. (|10|) and 
(|13p . The inset shows the ratio of the SQ and adiabatic expo- 
nents (diamond) and the n = 2 and 1 exponents (circle) from 
numerics, which agrees well with the theoretically expected 
value of 2. b.) Typical numerical results (solid lines) of the 
LE of the XXZ model are shown for a SQ from the XX point 
to J, — 0.5J and -0.3J, the dashed curve is the analytical 
expression using Eqs. ([9]) and c.) The scaling of the 

numerical data expected from Eq. H12p for short times for a 
SQ is visualized from = to O.IJ to 0.6J with O.IJ steps 
from top to bottom in arbitrary units. 



covering 1/2 < K < oo, for adiabatic ramps and SQ 
from = to a finite Jz . We performed the simulations 
using a combination of MPS 12 21 based infinite density 
matrix renormalization group |23l425l | to find the ground 
state variationally, and the infinite time evolving block 
decimation [2^ algorithms to simulate the quench. The 
LE is calculated by finding the dominant eigenvalue of 
a " generalized" transfer matrix which we obtain by con- 
tracting the tensors representing the two states. 

The factor L/na in the exponent in Eqs. PH)) and 
contains the unknown short distance cutoff. However, 
its value can be fixed by calculating the fidelity suscep- 
tibility, x/j around the non-interacting XX point of the 
Heisenberg model, in which case L/2Tra sa iVx/7r^, where 
N is the number of lattice sites and Xf ~ 0.0195[27|. 
Using then the Bethe Ansatz result [15[, K = 7r/2[7r — 
arccos( Jz/ J)], we find an excellent agreement with the 
numerical data with no fitting parameter (see Fig. [5^). 
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This excellent agreement is somewhat surprising since, 
as expected, the non-universal transient signals clearly 
differ in the LL approach and the numerics, and also, 
because the LL description completely neglects (asymp- 
totically irrelevant) back scattering processes, contained 
in the the lattice calculations. 

Slight deviations show up only at the end points of the 
critical region of the XXZ Heisenberg model in Fig. [21 
Close to the Jz = J point, the previously mentioned back 
scattering term, driving the Kosterlitz-Thouless phase 
transition causes a slight disagreement. Upon approach- 
ing the ferromagnetic critical point at Jr = — J, on the 
other hand, the validity of bosonization shrinks to very 
small energies, and the high energy modes, absent in the 
present considerations, also influence the overlap. 

Detection: We propose to measure the LE of the LL in 
a cold atomic setting, where a flux qubit 28, 2^ is used 
to control the interaction between a the atoms. Though 
here we focus on a one dimensional LL, the proposed 
setup also opens up the possibility to study the LE in 
higher dimensional interacting systems. The flux qubit 
consists of a Joscphson junction circuit, and is governed 



by the Hamiltonian H„ 



ibit 



ecTz + Aax, with A the 



tunneling between the two eigenstates of CTz, | O) and 
I O), carrying oppositely circulating persistent currents, 
±/, and e the energy splitting. In addition to the ex- 
ternal flux, <I>oxt, the states | O) and | O) generate an 
additional flux, Ideally, tunneling between them is 
suppressed. The one-dimensional quantum gas is pre- 
pared in a two-dimensional optical lattice potential (soj 
at the chip surface, hosting the flux qubit. A few mi- 
crons thick, high reflectivity (>99.9%) dielectric coating 
of the chip (not shown in Fig. [1]) can shield the super- 
conducting device from the laser power (mW) incident 
during the measurement time (ms). The quantum gas is 
positioned at the flux qubit (see Fig. [1]) by positioning 
the laser beams (optical tweezers), such that the total 
magnetic field of the state | O) of flux $oxt -I- be at 
a Feshbach resonance, while the field of the state | O) 
of flux ^ext ~ ^/j be further away from the resonance. 
We estimated the qubit switching-induced magnetic field 
difference for an elongated rectangular flux qubit, with 
parallel sides comparable to the length of a typical cold 
atomic tube (~ 10 /im). Assuming a persistent current 
of / = 2/iA and a separation oi d ~ 2R = 2fim between 
the two lines of the qubit, we obtain a field difference 
SBf 16 mG. Although relatively small, this field is 
comparable to the width AB = 15 mG of some narrow 
Feshbach resonances used to realize a LL in *^Rb sys- 
tems 



311. 



In this setup, one could use rf spectroscopy and mea- 
sure the absorption spectrum of the qubit [S^ in the pres- 
ence and in the absence of the trapped gas. Similar to 
the X-ray edge singularity problem, this absorption sig- 
nal is just proportional to the Fourier transform of the 
LE. In this case, the qubit can be positioned far away 



from the degeneracy point, at |e| ^ A, and the rf/mw 
spectroscopy can be referenced to high stability quartz 
oscillators {6f / f ^ 10^^'^) to resolve the fine structure 
at the absorption edge. 

Alternatively, the LE can be measured using Ramsey 
intcrfcrometry]^ Q: initializing the qubit in the | O) 
state with weak interactions to the cold atoms, yields a 
wavefunction | O) ® at i = 0. By applying a 7r/2 rf 
pulse (33! , a superposition of the two qubit states is pro- 
duced (I O) -f I 0))/\/2(8) l^'o), yielding distinct, qubit 
state dependent time evolution for |4'o). After time t, a 
second 7r/2 pulse and the measurement of the qubit cur- 
rent (/) ~ ((T2) is performed, giving a signal proportional 
to C{t). The timescalc separating Eqs. ((T^ and ([T^ is 
estimated as a/v 40 /xs from a typical Fermi tempera- 
ture of 1 /iK or trapping frequency of 25 kHz 3l|, and 
can be made even shorter by executing large quenches. 
Given the rapid progress of technology of flux qubits sj], 
these are already accessible with current coherence times 
and will be even more so in the future. 

Note added. After accomplishing this work, we became 
aware of related proposals to measure the work statistics 
in hybrid architectures (ssi [sij . 
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SUPPLEMENTARY MATERIAL FOR 
"LOSCHMIDT ECHO AND THE MANY-BODY 
ORTHOGONALITY CATASTROPHE IN A 
QUBIT-COUPLED LUTTINGER LIQUID" 



DERIVATION OF THE TIME EVOLUTION 
OPERATOR 



where the operators 



bp. 



(S2) 



K+{q) = b+bt^, K-{q) = b,b^, (S3) 

are the generators of a SU(1,1) Lie algebra, satisfying 
[K+iq),K_{q)] = ~2Ko{q), [Ko{q), K±{q)] = ±K±(q), 
and the operators for distinct g's commute with each 
other. 

Using these and a faithful matrix representation of the 
generators of SU(1,1) [l^, the coefficients in Eq. ([ST|) of 
a given mode are determined from substituting this to 
the Schrodinger equation of the time evolution operator, 
iU{t) = HU{t), or rather iU{t)U*{t) = H, where U{t) = 
dtUit). Using the SU(1,1) commutation relations, we 
derive the identities 

exp(aii:+)A'o = (i^o - aK+) exp(aii:+), (S4) 
exp(feA'o)A'_ = exp(6A'o - b), (S5) 
exp(aA + )A:_ = {K_ + a^K+ - 2aKo) exp(aA'+). (S6) 

Using these, we finally get 

^C+{q,t) + C+{q,t)Co{q,t)- 

-exp[~Co{q,t)]Cl{q,t)C-{q,t)^tg{q,t), (S7) 

tC-iq,t)cxp[^Coiq,t)] ^ giq,t), (S8) 

-Co{q, t) + 2C+{q, t) exp[-Co(g, t)]C. {q, t) = 2lLo{q, t), 

(S9) 

(pg{t) = uj{q,t), 



with initial conditions C+(g,0) = Co{q,0) = C ( 
V3,(0) = 0. These can be solved following Ref. [31 



C+{q.t) 



u*{q,ty 



C-iq,t) 



u*{q,t) 



(SIO) 

?,o) = 

as 



(Sll) 



C^{q,t) = -2\n{u*{q,t)), ^,{t) ^ J io{q,t')dt', (S12) 



where the time dependent Bogoliubov coefficients stem 
fromisi 



idt 



Ug{t) 









w(g,t) g{q,t) 
-9{q,t) -uj{q,t) 



Uq{t) 



with initial condition 



"«,(0)' 




■ 1 " 










(S13) 



(S14) 



The time evolution operator of a given mode is 

Ug{t) = exp[C+(q,t)A%(g)]exp[Co((?,i)i^o('?)]x 

xexp[C_(g,t)A'_(g)+z(p,(t)], (SI) 



BOGOLIUBOV COEFFICIENTS AFTER A 
BANG-BANG PROTOCOL 

Using the Bogofiubov coefficients after a SQ (i.e. Uq{t) 
and Vq{t)), the multiple quench scenario, shown in Fig. 



6 



can be analyzed. The time dependent Bogoliubov where 



coefficients at time t after the nth quench, u'q'\t) and 



(t) are determined from 



f/"(i) 



U{t) 



-iuJi {q)t 









1 


Uq{t) 

Vgit) U*Jt) 



(S15) 

(S16) 
(S17) 



where that last interaction change occurs for an nth or- 
der protocol at (2n — l)t as shown in Fig. IS1[ the time 
evolution afterwards only contributes with an irrelevant 
phase factor to the overlap, but it is important to retain 
it to check the convergence of the numerics. 
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FIG. SI. (Color online) The bang-bang protocol is visualized 
up to 71 = 3. The arrows indicate the time instant when the 
overlap is calculated for an nth order protocol. 



This leads to 



k")(t)p = i + 





2 







J2 T {s,{t) + ayg{t)f 



.<7 = ±1 



(S18) 



y,it) = - 1, s,it) = Re (u,(t)e— '(^)*) , (S19) 

Ug{t) = cos{ujf{q)t) - ^sin{ujf{q)t) + ^ 



(S20) 



and \vq{t)\'^ = \ug{t)\'^ - 1. By plugging Eq. §Wl to 
Eq. (9) in the main text, expanding the logarithm, then 
performing the momentum integral in the t ^ a/vf limit, 
and finally rcsumming of the resulting scries gives 

Cni2nt) ^ Q^it) ^ Cll (S21) 



ADDITIONAL PROPERTIES OF THE 
LOSCHMIDT ECHO 

From Eqs. (10) and (13) in the main text, the LE is 
found to be symmetric with respect to the initial and final 
LL parameters in the steady state {t — oo) for adiabatic 
and sudden quenches as {Ki,Kf) i — > (Kf,Ki) < — > 
{1/Ki, 1/Kf). Since these hold true in the two extreme 
cases, it is plausible to assume that they remain valid 
for any smooth, monotonous protocol. Throughout the 
calculations, we have also assumed L > vt, otherwise 
revivals show up with a period L/2v. 

Our adiabatic and SQ overlaps parallel closely to An- 
derson's orthogonality catastrophe and the related x-ray 
edge singularity, respectively. However, while the overlap 
decays exponentially with the system size in our case, the 
impurity problem yields a power law suppression of the 
overlap with the system size or time. 



